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MOTION OF A BALLISTIC MISSILE ANGULARLY MISALIGNED WITH 
THE FLIGHT PATH UPON ENTERING THE ATMOSPHERE AND ITS 
EFFECT UPON AERODYNAMIC HEATING, AERODYNAMIC 
LOADS, AND MISS DISTANCE 1 

By H. Julian Allen 


SUMMARY 


An analysis is given of the oscillating motion of a ballistic miss il e 
which upon entering the atmosphere is angularly misaligned with respect 
to the flight path. The history of the motion for some example missiles 
is discussed from the point of view of the effect of the motion on the 
aerodynamic heating and loading. The miss distance at the target due to 
misalignment and to small accidental trim angles Ib treated. The sta- 
bility problem is also discussed for the case where the missile is 
tumbling prior to atmospheric entry. 


INTRODUCTION 


It is characteristic of long-range rockets that, because of the low 
efficiency of the propulsion system, the weight at take-off is large com- 
pared to the final weight after fuel is expended. Typically, a saving 
of 1 pound in final weight can save of the order of 20 pounds of Initial 
weight and, as a result, strict attention must be given in the design of 
rockets to keep design safety factors to a minimum. Thus the magnitude 
of the factors which principally influence the final weight must be known 
with as great accuracy as possible. 

Two such factors are the aerodynamic load experienced by the vehicle 
as it descends through the atmosphere, which affects required structural 
weight, and the aerodynamic heating experienced in the descent, which 
affects required ceolant weight. Problems relating to the loading and 
heating of missiles during atmospheric entry have, of course, been given 
considerable attention, both from a general point of view (e.g., ref. 1) 
and in detail for specific designs. In the usual treatment of the problem, 
however, the rather idealized case has been treated wherein the vehicle 
enters the atmosphere unyawed or unpitched with respect to the flight 
path and without angular velocity. If the vehicle enters the atmosphere 
in a yawed or pitched attitude, it will, during its oscillatory approach 
to the earth, be subjected to lateral forces in addition to the longi- 
tudinal forces due to deceleration. Moreover, the distribution of 
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aerodynamic heating over the surface for the oscillating vehicle will 
differ from that for the vehicle if aligned with the flight path. Ihus 
a question arises as to what extent the structural weight and the weight 
of coolant might -he altered hy the fact that the rocket upon entering 
the atmosphere is angularly misaligned and has angular velocity. 

The analysis of reference 2 provides an excellent basis for such a 
study. However, the results of that analysis are in a form which is not 
convenient for demonstrating the relative importance of the several fac- 
tors which are of principal interest to the loading and heating problems. 
It is the purpose of this paper to re-examine the motion analysis of 
reference 2 using some simplifications which were employed in reference 1 
in order to indicate more clearly the salient features of the motion 
problem and its effect, in turn, on loading and heating. 2 


BASIC ASSUMPTIONS 


In the analysis to follow it will be assumed that the missile warhead 
which enters the atmosphere is rotationally symmetric so the misalignment 
angle may be considered as yaw or pitch or any vector combination thereof, 
and that the fineness ratio is sufficiently low and the Mach number is 
sufficiently high that the pressure distribution is independent of Mach 
number for the Mach number range of importance (see ref. 4). Thus, the 
rates of change of the respective aerodynamic force and moment coeffi- 
cients with a, d, and q are considered to be constants. The basic 
assumption of the analysis of reference 2 is retained in the analysis to 
follow; namely, that the angular oscillations are small so that the sine 
of the angle of oscillation is the angle of oscillation in radians and 
the cosine is unity, and the drag coefficient is sensibly the same as it 
would be for the nonoscillating missile. In addition, the assumptions 
employed in reference 1 are also employed herein; namely, that the accel- 
eration of gravity is constant with altitude, the flight path through the 
atmosphere, is essentially a straight line, and the variation of air density 
is the exponential function 


P 



( 1 ) 


wherein p Q and p are constants and y is the altitude measured from 
sea level. 

A complete list of symbols is given in Appendix A. 

2 It should be noted that an analysis has been made in reference 3 
which, although principally aimed at study of other features of the sta- 
bility problem, employs a basic approach similar to that of the present 
report. 
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ANALYSIS 


If the angular displacements are small, then the differential equation 
of angular motion with time as the independent variable may be written 


fit + f i(t) + fa(t)a = 0 

dt dt 


( 2 ) 



wherein the time-dependent coefficients are given by 3 


fi(t) 


C^pVA (C^ + C m4 )pVAZ 2 
2m 21 


(3a) 



’C^p V 2 AZ 
4lm ’ 21 


(3b) 


where the angle of attack is as indicated in the sketch. 

s This formulation is equivalent to equations (12) and (13) in ref- 
erence 2 except that in reference 2 the value of has been tacitly 
assumed to be zero, which is a justifiable assumption at high Mach num- 
bers. The quantity is retained in equation (3a) to be consistent 
with NACA standard nomenclature (e.g., see ref, j?). 
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It is convenient , now, to rewrite equation (2) with altitude, y, 
rather t han time, t, as the independent variable. To this end, it is 
noted that for the straight flight path assumed (0 = % = constant) 



-V sin 0 E 


(*a) 


d 2 y 

dt 2 


52 al n e E - v 52 

dy dt dy 




Thus 


do, = dn dy 
dt dy dt 


-V sin 



(5a) 


d 2 q _ d 2 a f dy^j 
dt 2 dy 2 


dq d 2 y _ 
dt 2 ~ 


V' 


2 sin 2 0- 


E 


d 2 q 

dy 2 


V Bin 2 e E ^ 
dy ^ dy 


(5b) 


so that equation (2) becomes 


d 2 q 

dy 2 


+ f*3(y) ^ + f*(y)« 


0 


( 6 ) 


wherein the altitude-dependent coefficients are 


fs (y) 


dV/dy 

V 


C LqP A 


(Cmq + Cjh^) pAZ 2 
21 sin 0 E 


(7a) 


2m sin 0 E 


+ 
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C^ACdp/dy) C^ApCdV/dy) C^C^A 2 ! 2 c^pAl 
2m sin 0 E 2mV sin 0 E 4lm sin 2 © E 21 sin 2 ^ 


(7b) 


It is shown in reference 1 "that by use of the same basic assumptions 
as have been made in this report, the velocity may be expressed as 


V = V E e 


ko -py 

-5- e 


(8a) 


where 


^o - 


c DPq A 
pm. sin © E 


(8b) 


Hence 


dV /dy _ Pko -py 
V 2 6 


8(c) 


and from equation (l) 


_ -py 

dy " Pp o e 


(8d) 


Thus equations (7a) and (7b) become, upon setting the square of the radius 
of gyration as cr 2 = i/m. 


_ , , PqA r ( c mq + Cm^.)^ __ 

fs (y) = On - Cr + — — - — e py 

2m sin 0g L U ^ a 2 


(9a) 


f4r(y) =■ 


PqA c^i 

_ 


5TS^L"~ +C V sln M 


e Py t Po 2 A 2 


4m 2 sin 2 0 E L 




Cm q C La Z 


2-, 


,- 2 Py 


(9b) 
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If we set 

then equation (6) may be written 

-Y 


Y = |3y 


d 2 a 

dY 2 


+ 2kj_e " ■— + + ^3 e = 0 


(10a) 


(10b) 


wherein the constants are 




kn = 


k 2 " 


k3 - 


Po A 


C D “ C I tt + 


' Offiah 


4pm sin 9^ 

Po^ _ 

2p 2 m sin 2 0g V cr 2 

Po' 


(Cmq + c i%)(o) _ 


+ C^p sin 0g 


2 A 2 


4p 2 m 2 sin 2 0g 


“ ( V 3d - J 


J 


(10c) 


In order to find a solution for the differential equation (10b) let 

-Y 


k^e 


so that 


a = rje 


to _ k ie " Y dq -Y^e"* 

ar e dY kie e n 


(ii) 


and 


da 

dY 2 


k ie - T ^ 
dY 2 


sa + fi 

dY V 


-Y -Y- 

-Y k^e . 2 -zY k x e 
k z e e + k x 2 e e 


) 


Substitution in equation (10b) then yields 
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~ ~ 2 + (^-2 + ^ + (^3 ” k x 2 )e 2 

dY L 


T] = 0 


( 12 ) 


The rigorous solution of this equation is not necessary since the 


term 


10 


(ka “ k x 2 )e“ 2Y 


is, for the type of missiles to he considered later, small compared to 

,-s 

|0 

(k 2 + k x )e~ Y 


particularly at the higher altitudes which are of principal concern. If 
the square term is omitted, the solution of equation ( 12 ) is known (see, 
e.g. , ref. 6 ) 





+ k x e 



+ C 2 Y, 


0(2^ 


+ k 3 



(13) 


where C x and C 2 are constants of integration and, in accordance with 
Watson’s notation (ref. J) , the functions J Q and Y 0 are the zero order 
Bessel functions of the first and second kind, respectively. 

Combining equations ( 10 a) , (ll) , and (13) gives finally 


k x e 

a = e 


■Py 


'iJq^2 \f 


k 2 + k x e 


*) 


+ c 2 Y ( 


o^2 \f 


k^ + k x e 


3y- 

2 


) 


(14) 


If 

missile 


it is specified that on entering the atmosphere (y -> 00) the 
is misaligned by the angle og but has no angular velocity, then 




C x - ctg 

c 2 = 0 


(15a) 
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and equation (l4) can be written 

— = e klS Jd^2^kg + k x e "^) (15b) 


In the cases of usual interest the quantity k 2 is very much larger than 
kj., as will be shown later, so that one may use the approximation 





(15c) 


For large values of the argument, the approximation 


a 



(154) 


can be used, and the maximum or envelope value is thus 



k x e 


■py 




s[k2 


py 

2 


(I5e) 


In Appendix B are derived expressions for the angular velocity and 
angular acceleration as well as their maximum values during any cycle of 
oscillation. The approximate maximum angular velocity per cycle is (see 
eq. (B5c)) 


"d(a/ag) 

dt 



BY-gSin 9ge 




( 15 *) 


while the maximum angular acceleration per cycle is 
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rd 2 (a/^)f k 2 ‘ s , „ -2|z (ki-koJe'^ 

rs?-J = /; pv ® w v e 


(15b) 


From equation (l5d) it can be found that one cycle of oscillation 
takes place during the altitude change 


2rt 

. 3y 

3 n/ k 2 e 2 

The frequency of oscillation is thus 


_ pz 

P^/k^ e 2 V sin e E 

to c — , — — 

2rt 


and using the velocity from equation (8a) one has 


Ci] = 


pv EN /k^sin e E _J^ e "Py_4i 


2it 


(16) 


The frequency is maximum at the altitude 


yi = - in ko 


(17a) 


which, as shown in reference 1, is the altitude at which the deceleration 
due to drag is maximum and the velocity is 


_ 1 

V x - e 2 V E E 0.6l V E 


(17b) 


The maximum frequency is 


“max = “1 = 


^ /'$ v E e 2sln GgN 

tK s — ) 


(17c) 
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and the corresponding amplitude is from equation (I5e) 



(174) 


For later use it is of importance to note that in the analysis of 
reference 1 it is shown that for turbulent flow the altitude for which 
the average heat rtransf er rate .is a maximum is 


y 2 = 


l 

3 



when 


_ i 

V 2 = e" 3 V E = 0.72 V E 


(18a) 


(18b) 


The corresponding frequency is 


Wg = 



(18c) 


and the corresponding amplitude is 



On the other hand, for laminar flow (which, it is expected, should 
be applicable at least for the stagnation point on the missile nose), the 
altitude for' maximum heating rate is then 


73 =--£ Jn(3ko) 
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when 


and 


and 


_ i 

V 3 = e" 6 V E = 0.85 V E 



( 19 b) 


( 19 c) 


(I9d) 


From a loads point of view it is of value to note that from equa- 
tion ( 15 g) the altitude for which the maximum angular acceleration 
occurs is 


=| l*\ 


0 (ko - *i) 


(20a) 


and that the maximum acceleration is 


'd a (q./a E ) 
L dt 2 


•nax 


= = _i_ r 

L dt 2 J* Jx u 


3*2 


M*o - *i) J 


4 P 2 V E 2 sin 2 e E 


(20b) 


Also the normal force experienced is 


N = \ C Na aV 2 pA 


(21a) 


and substitution from equations (l), (8a), and (l5e) gives, for the maximum 
normal force per cycle 


N 


* _ 


c Hq a EPo v E gA e -f py e (fc 1 - ko)e" Py 
2 \j Tt *J ~ *2 


(21b) 
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This force is a maximum at y 4 and has the value 


N, 


max 


= H 4 *- 


C Ng, a EPo V E A 
2 si n /"^2 


{[«ano]‘ *> 


(21c) 


Finally, it is of interest to determine the order of magnitude of 
the drift due to Og from the course the missile would have had if aligned 
with the flight path. If n is defined as the distance normal to the 
straight line trajectory (that the missile would have if . ctg were zero), 
the lateral acceleration is 


dfn 

dt 2 


-V sin 9- 


E 


<3y 


/W\ = L 
\dt/ m 


C I<x aV 2 pA 


2m 


( 22 ) 


Use of the relations of equations (l), (8a), and (15c) gives the lateral 
velocity at altitude y 


dn C hx a E V EPo A r Y 


dt 2m sin 6- 


E 


l 






Py- 


') 


<3^1 (23a) 


With the substitution 


?! - 2^/^ e 


.IJl 


then this may he written 




P7 


dn = c L a a E V EPo A 
^k 2 0m sin 9- 


E 


I 


eN 4*2 J 1 5 x Jo(5 1 )dg 1 


(23b) 


The lateral displacement at y ^ 0 is, by similar substitution, 
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n = 


2 sf^z /^ 0_'\- 2 
C Ict a EPo A £ e v sk 2/ r 


2 k 2 p 2 m ein 2 0 g 




k . 


^JoCSjcL^ (23c) 


Defining 


F(ko,kx,k 2 ) 


2«/k^ 

s i_ r e^£_ r 
k 2 4 5 Vb 



£ 1 kl ~ a j £ 2 

e \ ^ 2 / 1 £ i J 0 (£ i )d£ 1 d£ 


(24a) 


then 


Clo^Po^ „„ , , , 

n = F(ko,kx,k 2 ) 


2 p 2 m sin 2 0 - 


(24b) 


E 


The integration of equation (24a) may be performed in the special caBe for 
which k x = ko = 0 , for then 

2*/^ 


*0,0,1*) 


•/kg* 


■5 / 


Jx(£)d£ = 


1~ J~o(2 n/ kg) ^ 

kg k 2 


For values of k Q and k x other than zero, an analytic solution of equa- 
tion (24a) is not known. However, Dr. William Mersman has dete rmi ned 
values of the function by numerical integration using an IBM 65 0 type 
digital computer. The computational procedure and results are given in 
Appendix C. Over the range of variables of interest in this paper, 

0 < ko < 10 , -80 < k x < 0 , 5 xl 0 4 < k 2 < 8 xl 0 5 , the integration indicates 
that 

F(ko,kx,k 2 ) = l/k 2 

within 0.2 percent (see -Appendix C, eq. (C20) ) . 
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If the angle og is considered as a yaw angle then the miss distance 
is the lateral or "deflection target error" 



C I <x (X EPo A 
2p s m sin a 0 E 


F(ko,k!,k 2 ) 


(25a) 


while if the angle cog is considered as a pitch error , then the miss 
distance is the longitudinal or "range target error" 



C lc L a EPo A 
2p 2 m sin 3 0 E 


^(kok^ks) 


(25b) 


Thus, in the general case, the area of miss is elliptical with the 
major axis in the range direction (except when the missile descent is 
vertical when the miss area is circular) . 

It is useful, for comparison purposes, to find the miss distance which, 
results from an accidental trim angle, Orp. The differential equation (22) 
becomes in this case 


d 2 n 

dt 2 



c I < x a T? vg P A 

2m 


(26a) 


so that 


an = _ C WPq V E A f 7 e - ^ e^i ^ 

at 2m sin 0g ^ 1 


or 


dn 

dt 


c I a , a TPo v E A / 
koPm sin \ 



and, in turn. 


= C I<i co TPo A r°f 1 

koPm sin 2 0 E ^ V 



(26b) 
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which, can "be shown to give 


n 


Cl^PpA 

P 2 m sin 2 Sg 


Fi(ko) 


( 26 c) 


wherein 


Fi(ko) = 


ho 


Ei^) -ln(^) - 0.577216 


( 26 d) 


and ET(k Q /2) is the exponential integral for which tabulated values are 
given in reference 8. Values of P 1 (k 0 ) have been computed for practical 
values of ko and are presented in table I. 

If the trim angle is a yaw angle then the miss distance is the 
deflection target error 


(^t)* ■ fftn Fl(ko) 


(27a) 


while if the trim angle is a pitch angle then the miss distance Is the 
range target error 



G i a a rPo A 
P 2 m sin 3 0-g 


Fi(ho) 


( 27 b) 


DISCUSSION 
Conical Warheads 


It is the purpose in the discussion to follow to examine the angular 
motion of typical ballistic missiles In the atmosphere to ascertain the 
importance of this motion to the problems of aerodynamic loading and 
heating and miss distance. Conical shapes have been chosen for this study 
because the calculated stability derivatives are available. 

Tobak and Wehrend (ref. 9) have calculated the stability derivatives 
for cones of half -angle, 6. Although they give results which are applicable 
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from low supersonic speeds (exceeding the Mach number for shock detachment) 
to hypersonic speeds, the concern of this paper will be only with the 
hypersonic or Newtonian solutions. For arbitrary distance from the cone 
apex to the center of gravity, 2 C g, reference 9 gives the pertinent 
derivatives, using the symbols of this report, as 




Cjy = 2 cos 2 & 


C x = C D = 2 sin 2 5 

= C Wa - C D = 2(cos 2 S - ein 2 5) 



^ + 2 cos 2 8^-y^^ 



J 


(28a) 


In addition it should be noted that for cones the center of volume is at 
32/4 and the square of the ratio of cone length to radius of gyration for 
arbitrary center-of -gravity position is (considering the body to be 
uniformly solid) 



80 


12 tan 2 5 +3+80 



(28b) 


Validity of the Analysis 


Before examining the stability, loading, heating, and miss-distance 
problems it is necessary to determine whether the previous analysis is 
valid for the conical shapes to be considered. It was noted that the 
solution given by equation ( 15 b) applies only in those cases for which in 
the differential equation (12) the term 
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(ks - k 1 2 )e" aY (29a) 

can be neglected in comparison with 

(k 2 - k 1 )e" Y ( 29 b) 


In Appendix D it is shown that the values of the factor (29a) are, 
for practical cases , always very small compared to the values of the 
factor ( 29 b) and, in addition, that ks is, for practical cases, always 
much larger than k x so that the solution of equation ( 15 c) is valid. 


Stability 


In the equation (15c) it is clear that the missile is statically 
stable if k 2 is positive and certainly dynamically stable if k x is zero 
or less than zero. If k x is positive, then it is possible for the oscil- 
lation amplitude to increase with decreasing altitude. Thi s may be con- 
veniently shown from the approximate equation (I 5 e) for the envelope value 
of a. Writing this equation in the form 



(30) 


we find that the derivative of this function with respect to y is zero 
when 



or 


y = | In^ki) 


(3D 


Thus if positive ki is even as large as 0.25, the analysis indicates 
that the incipience of divergence occurs (at sea level). For larger posi- 
tive values serious divergence at the lower altitudes would be anticipated. 
In order to obtain a better grasp of the nature of such motions it is 
instructive to examine the indicated behavior of a missile with some arbi- 
trarily assumed (but practicably realizable) static stability for several 
values of k x . For this purpose let it be supposed that k 2 is 10 4 and 
k x has, in turn, the value -10, 0, and +10. The equation (l5e) 
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then yields the envelope curve of angular history shown in figure 1. It 
is indicated that positive values of ki could promote serious divergence 
of the amplitude of .angular oscillation. Friedrich and Dore (ref. 2) dis- 
cussed the possibility of such divergence and noted that these adverse 
effects could occur if the missile underwent large reductions in speed 
during the descent. However/ no consideration was given in their report 
to the importance of the damping terms. In other words what they con- 
sidered, in the language of this report, was that the coefficient k x (as 
given by eq. (10c) ) was overwhelmingly influenced by the drag coefficient. 
These results of their analysis thus imply that high drag shapes are unsat- 
isfactory for ballistic missile application in spite of their inherent 
advantages in the aerodynamic heating problem (see, e.g., ref. 1). 

The question of importance is, then, whether or not it is realistic 
to ignore the damping terms (Cj^) an< 2- (C m ^ + C^) due to plunging and 

rotation, respectively, in the determination of sign and magnitude of 
ki. To answer this question it is convenient to consider a simple conical 
shape of arbitrary cone half -angle, 5. 

To investigate the sign of ki for cones it is sufficient to 
evaluate the "dynamic stability" factor 


C D - + (C^ + 0^5 


in equation (10c) by use of the relations of equations (28a) and (28b). 

This has been done and the results are presented in figure 2 for several 
cent er-of -gravity positions. (Note that l C g/l = 3/4 is the center of 
volume and a most likely position far the center of gravity. ) From fig- 
ure 2 it is seen that ki must always be negative for conical (and pre- 
sumably for near conical) shapes. Thus the inference of reference 2 that 
high-drag (i.e., large 5) shapes are necessarily undesirable from the 
dynamic stability viewpoint is unjustified. Nevertheless, "it should be 
noted that while for the high-drag cones previously discussed the value 
of k x was _always negative, high-drag shapes which are blunted at the 
nose can, in fact, have positive values of kj.~ This leads to some concern 
since, for reasons of aerodynamic heating, considerable blunting of the 
nose may be desirable. However, it should be noted that the divergences 
indicated by the an alysis given previously are greater than what will 
actually occur for the following reason : The analysis employs the velocity- 

altitude relations from reference 1 and so does not include the effect of 
gravity on the velocity-altitude history. This assumption is admissible 
for the high-speed portions of the trajectory and for the portions wherein 
the deceleration is large compared to that of gravity. For low drag s ha pes 
the assumption is generally admissible over the whole altitude range but 
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for high drag shapes it becomes inadmissible at low altitudes. Thu s the 
actual trajectory is not a linear one but is steepened at the low alti- 
tudes and the actual velocity approaches the so-called "te rmina l velocity" 
(i.e., that speed for which the drag equals the weight) rather than zero 
velocity. Then, in actuality the term (dv/dy) /V in equation 7(a) 
approaches zero rather than the value indicated by equation 8(c). Hence, 
while at high altitudes kj_ is determined by the "dynamic stability" 
factor given above, at low altitudes for high drag shapes it would be 
determined more nearly by 


- C La + 



and thus the actual value of k x at low altitudes would be more negative 
than the analysis would indicate. Therefore, while divergence of the 
motion might occur for a high drag shape, it would not necessarily 
continue . 

Next it is in order to examine the sign of the factor which controls 
the static stability parameter 


k 2 = 


Po A 

2p 2 m sin 2 0g 



sin 9 



It is shown in Appendix D that the second term in the bracketed expression 
Is small compared to the first for practical cases so that it is only 
necessary to be assured that is negative to insure stability. In 



figure 3 the Newtonian value of 


C. 


(from eqs, ( 28 a) and ( 28 b)) is 
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plotted as a function of 5 for various values of l C g/l, and it is seen 
that when S is small, care must "be exercised to keep the center of gravity 
far forward. It should he noted (see ref. 9) that the center of pressure 
is independent of Mach number down to the Mach number of shock detachment. 
Thus the hypersonic requirement of l C g is also the Supersonic requirement. 

From the preceding discussion, it is apparent that at least from low 
supersonic to hypersonic speeds, positive static stability (i.e., positive 
k 2 ) can be obtained. Similarly, it can be shown that over the same speed 
range dynamic stability is assured (i.e., negative k x ) . Now it is impor- 
tant to determine the magnitude of the static and dynamic stability which 
can be provided. To this end the following digression in the discussion 
is in order. 

For long-range .ballistic missiles the aerodynamic heating problem 
must be the principal consideration in design. It has been shown (e.g., 
see ref. 1) that, generally, the aerodynamic heating problems are 
reduced when the value of ko is increased. On the other hand, if ko 
is too large then the speed of descent of the_ missile becomes low for too 
great a part of the final trajectory which increases the vulnerability of 
the missile. Thus some compromise is required and this compromise value 
of k Q tends to be larger the longer the range. A value that will be 
considered herein to be a reasonable one for a 3,000- to 5,000-mile range 
would be of the order of 5 to 20. 

Now, k x (see eqs. (8b) and (10c)) can be written in terms of ko in 
the form 


ki 



( 32 a) 


while for kg, if the 0^(3 sin 9 E part is neglected (see Appendix D), this 
parameter becomes 


jfo f ’if 

2 (31 sin e E vW V V 


( 32 b) 


Values of k^/ko and [31 sin 9gk 2 /ko are given in tables H and III. Since 
k x /k Q depends upon the location of the center of gravity and the cone 
angle while k 2 /ko depends upon these factors and, in addition, the length 
of the missile, it is necessary to consider some examples in order to 
determine the magnitudes likely to be realistic for k x and k 2 . Accord- 
ingly, let it be assumed that the missile weight is 3,000 pounds, that the 
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entrance angle of the trajectory, frg, is 30 °, and that the values of the 
atmospheric density relations are those of reference 1 (p Q = 0.0034 slugs 
per cubic foot, = 22,000 feet). For values of ko equal, in turn, 
to 5, 10, and 20 the base diameters, lengths, and volumes of the example 
missiles are those given in figures 4(a), 4(b), and 4(c), respectively, 
as functions of cone half-angle. For this analysis it is arbitrarily 
assumed that the maximum allowable missile length is 30 feet and the mini- 
mum allowable volume is 10 cubic feet (corresponding to a high missile 
density of 300 pounds per cubic foot) . Thus in figures 4 the curves extend 
only to the cone half-angles which correspond to these two limits (the 
small cone-angle limit corresponds to the maximum allowed length and the 
large cone-angle limit to the maximum allowed density) . Xn addition, it 
is arbitrarily specified that the center of gravity in each case is at a 
distance from the apex (leg) where the local diameter is 2-1/2 feet. The 
resulting ratios of l cg /l are shown in figure 4(d) . 

With these physical characteristics, the values of k x and k 2 are 
those of figures 5 and 6. It is seen that the dynamic stability is great- 
est for large values of the drag parameter but for small values of the cone 
angle. On the other hand, the static stability is generally greatest for 
large values of both the drag parameter and cone angle. A notable excep- 
tion to this trend of the static stability parameter is the sharp decrease 
of k 2 at the largest angle for the ko = 5 case. This sudden reduction 
results from the fact that the center- of -gravity position has rather 
closely approached the center of pressure. 

Consider, now, two extreme cases: first, the ko = 5 missile at maxi- 

mum allowed density and second, the ko = 20 missile at maximum allowed 
length. The former has least values for both stability parameters and 
therefore will oscillate with the largest amplitudes, while the latter has 
the largest dynamic stability parameter and has a rather high value for 
the static stability parameter and thus should be representative of the 
opposite extreme. The angular behavior with altitude for these two mis- 
siles has been calculated using equations ( 15 c) and (l5e) and is shown in 
figures (7a) and (7b). The high altitude oscillations of figure 7(a) are 
similar to those of (7b) but displaced downward, altitudewise, by about 
25,000 feet. This is an effect of the lower static stability parameter 
for the k Q = 5 missile. At the lower altitudes the ko = 20 missile 
oscillations decrease more rapidly by virtue of the larger dynamic stability 
parameter . 


Heating 


It was pointed out in the Introduction of this report that when the 
time rates of aerodynamic heating are largest, it is important that the 
oscillation amplitudes be small in order that additional coolant mass will 
not be required to protect the vehicle from excessive local heating. It 
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is to "be expected that, at least for not -too -small cone angles, maxi mum 
oscillation amplitudes of the order of a few degrees of arc should he 
permissible with no important adverse effects. 

To determine whether or not the oscillations will be important as 
regards aerodynamic heating it is again convenient to consider a particular 
example. The same 3 , 000-pound missiles are used in this study. Since the 
laminar heating rate always reaches a maximum at a higher altitude than 
does the maximum for turbulent heating, it follows that (a/ct<g)* at maximum 
heating will be greater in the laminar case . - The amplitude ratio at the 
altitude for maximum laminar heating rate (calculated using eq. (I9d)) jLs 
shown for the example missiles in figure 8. It is seen that these values 
are so low that no complications of the maximum heating rate problem due 
to initial angular misalignment with the flight path, Og, should be 
expected. In figure 9 both the angular amplitude ratio and the ratio of 
laminar heating rate to maximum laminar heating rate (Bee ref. 1) are 
plotted as a function of altitude for one particular example (ko = 5, 

S = 25°), which shows that while' the amplitude ratio is very small at the 
altitude for which maximum heating occurs, it may become sufficiently 
important at the higher altitudes where the heating . rate is still fairly 
high to require consideration in design. 


Loads 


To show the degree to which lateral loads due to ag are important 
it is again useful- to consider the example missiles considered earlier. 

To evaluate the maximum normal force using equation (21c) it is necessary 
to specify the entrance speed, Vg, and it will be assumed, for the exam- 
ples, that the speed is 20,000 feet per second. The maximum normal forces 
in terms of missile weight per degree angle misalignment at atmospheric 
entrance for the examples are shown in figure 10. While the normal forces 
are increased for. the longer missiles (due to .increased surface area), it 
does not appear that they could be too serious in a practical case. A 
20° value for oog only promotes a 3g normal acceleration for the long 
missiles which is small compared to the deceleration due to drag which 
(from the analysis of ref. 1) is 51g. Moreover, the maximum normal loads 
are not additive to the maximum drag loads since, as seen in figure 11, 
they occur at different altitudes. 


Miss Distance 


Before considering the actual magnitudes of the miss distances due 
to ag or ccrp it is well to discuss the accuracy of the analysis of miss 
distance given previously. In the analysis" it is assumed that the velocity 
at, all points of the trajectory is given by the exponential expression of 
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equation (8a.) and this expression was, in turn, obtained by neglecting the 
effect of gravity. As noted earlier the neglect of the effect of gravity 
is uni mportant in the evaluation of the velocity-altitude history except 
when the velocity is low and, simultaneously, the deceleration due to 
drag becomes comparable to the acceleration of gravity. Then the velocity 
given by the analysis falls below that which would act ually occur. In 
the analysis of miss distance it should be clear that the miss distance 
increases rapidly as the velocity decreases, and, hence, the miss distances 
given by the analysis are in error by an amount which increases rapidly 
with increasing when the decelerations near y = 0 become of the order 

of the acceleration of gravity. To assure that the deceleration at sea 
level is not less than lg, it is required (from ref. 1) that 



< Ss. 

sin 0^1 


For example, for an entrance speed of 20,000 feet per second and for 
©E = 30 °, the sea level deceleration reaches lg for ko of about 7 (see 
ref. l) ; hence, the values of F(ko,k!,k 2 ) (see Appendix C) and Fi(ko) 
from table I should not be used, under these conditions, for values of 
ko in excess of about 7, particularly as ko greatly exceeds this value. 
Thus in the calculations to follow the miss distances for k Q = 20 are not 
included and, in addition, the reader must note that even for ko = 10 the 
estimated miss distances exceed the actual ones. 

The range target errors per degree angle misalignment at atmospheric 
entrance for the example missiles (Vg = 20,000 ft/sec, 6g = 30 °) are shown 
in figure 12. The deflection target error, not shown, is simply one-half 
the range value. It is seen that, as with the normal force, the miss 
distance is greatest for the smallest cone angles. However, the miss dis- 
tance is trivial since even for an ccg of say 20° the range target error 
is but about 20 feet in the worst case. 

A serious problem is the miss distance which will result from a trim 
angle even slightly different than zero. In figure 13 is shown the range 
target error per degree of trim angle for the ko of 5 and 10. It is 
seen that a trim angle of as little as 0.1° can cause a range error of 
many miles. As the cone angle increases then for a given value of ko, 
the miss distance diminishes until when the cone half-angle is 45° the 
miss distance is zero since the lift-curve slope is then zero (see 
eq. (27a)). Except in this special case, however, the miss distance due 
to even a slight trim angle is very important. One method for reducing 
this miss distance would be to spin the missile about its axis so that 
it would follow a corkscrew path during descent. This solution introduces 
another difficulty, however, in that care may have to be taken to keep the 
spin rate from approaching the pitching (or yawing) frequency else tumbling 
may occur if the missile is not identical as regards aerodynamic and 
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inertial characteristics about any radial axis (see ref. 5 or 10). 
Unfortunately, the pitching frequency varies from zero to the maximum 
value given by equation (l6). The maximum frequencies for the example 
missiles are shown in figure 14. One obvious way to avoid tumbling 
resulting from "roll coupling" would be to spin the missile at a rate 
which exceeds, by a good margin, the maximum shown in figure 14. If the 
spin rates were the maximum pitch rates, the rim speeds at the base (i.e., 
at maximum diameter) would be those shown in figure 15. Since the required 
spin rate would have to materially exceed this rate, it is clear that a 
serious stress problem due to centrifugal loading might result (especially 
for small cone angles). In consequence, a better solution might be to 
spin the missile at a rate which is always less than the value of the 
pitch frequency at any given altitude, but the difficulty, then, would 
be one of assuring that the spin rate could hot accidentally approach the 
pitch frequency. In the preceding discussion it has been tacitly assumed 
that the missile is not identical as regards the aerodynamic and inertial 
characteristics about any radial axis so that the spin rate must not at 
any time match the pitch rate. Since the pitch rate changes rapidly with 
time (particularly at the higher altitudes), it is probably not a justi- 
fiable requirement that the spin rate not ever be the pitch rate, particu- 
larly since the asymmetries which exist may be trivial. Some further 
consideration must clearly be given this problem. 


Effect of Initial Tumbling 


In the discussion to this point it has been assumed that the missile 
enters the atmosphere misaligned by an arbitrary but fixed angle with 
respect to the flight path. When the missile is actually tumbling before 
entering the atmosphere, then the analysis given previously cannot be used 
since the equation of motion is restricted to small-angle considerations. 
(That the analysis is inadmissible is reflected in the fact that if, in 


equation (13 ) , C 2 


has a value which is not zero, then 


C 2 Y q (2 s/ks+kx e" 2 ^) 


becomes infinite if y is infinite.) In spite of this deficiency, some 
general remarks can be made about the effect of initial tumbling. It is 
clear at the outset that the missile must have but one possible trim 
attitude if initial tumbling occurs. If not, it could descend at some 
attitude for which no adequate protection for aerodynamic heating and 
loading had been provided. Furthermore it is a requirement that it must 
be righted to about the correct attitude at an altitude which is suffi- 
ciently high th'at the angular motions will become small by the time the 
heating and loading are intense. One obvious way which might assure that 
these conditions will be met would be to dispose of the empty fuel and 
oxidizer tankage only after the missile has entered the atmosphere 
sufficiently far to adequately correct the attitude. 
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CONCLUSIONS 


Prom an analysis of the motion of a 'ballistic missile initially mis- 
aligned with respect to the flight path prior to the entry into the 
atmosphere, it is concluded that it is possible to: 

1. Provide a continuously damped oscillation history with descent 
through the atmosphere. 

2. Keep the oscillations to a small amplitude when at altitudes for 
which aerodynamic heating and loading are severe. 

3. Prevent excessively large loads due to the oscillating motion. 

Moreover, while the miss distance at the target due to the initial 
misalignment angle is trivial, the error that can occur due to the trim 
angle being even slightly different than zero can be very large and its 
effect must be minimized 'in some manner. 

Since tumbling may occur prior to entry into the atmosphere the mis- 
sile must have only one trim attitude and must be brought near this attitude 
before the missile has progressed too far down through the atmosphere. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif . , 


for Aeronautics 
June 15, 1956 
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APPENDIX A 
SYMBOLS 


constant (See Appendix C, eqs. (C2).) 

reference area for coefficient evaluation (base area 
for cones) 

constant (See Appendix C, eqs. (C2).) 
constant (See Appendix C, eqs. (C2).) 
constants of integration 
drag coefficient — 

rate of change of lift coefficient with angle of 


attack, f -s^- 


a 0 


rate of change of moment coefficient with angle of 


attack, ( — 


a -> 0 


change of moment coefficient with time rate of change 

8c m 

of angle of attack, 


\ 


A 0 


rate of change of moment coefficient with angular 
8 Cm 

velocity, 


v) 


q. 0 


rate of change of normal- force coefficient with angle 
3Cjj 

of attack, — 

acta, 0 


axial-force coefficient 
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d 

e 

fa(y) ^4(y) 

Fi(k 0 ) 

F(ko,k!,k 2 ) 

S 

h 

I 

I 0,l,2 J1 , . 

J 0 ( ) 

J l,2,3,. . . ,r,s( ) 

ko 

ki 

k 2 

^3 

l 

l os 

L 

m 

n 


N 


diameter of "body base 
Naperian base 
functions of time 
functions of altitude 

function used in evaluating of "miss distance" due 
to cCfp (See table I.) 

function used in evaluation of "miss distance" due 
to ctg (See Appendix C.) 

acceleration due to gravity 

integer (See Appendix C.) 

mass moment of inertia 

integrals (See Appendix C.) 

Bessel function of the first kind of zero order 

Bessel function of the first kind of order 1,2,3* • • . >r,s 

the "drag" parameter (See eq.. (8b).) 

the "dynamic stability" parameter (See eq.. (10c).) 

the "static stability" parameter (See eq.. (10c).) 

the "cross-products" parameter (See eq.. (10c).) 

body length and reference length for moment coefficient 
evaluation 

distance from body bow to center of gravity 
cross-wind force 
missile mass 

distance normal to the trajectory the missile would 
have if it were angularly aligned with flight path 
and without angular 1 velocity 

normal force (force perpendicular to the axis of 
revolution) 
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P i,2,3,...,r,s( ) 

q. 

r 

s 

S 0,l,2,.. . ,T,b( ^ 
t 

V 

V E 

x 


functions (See Appendix C.) 
angular velocity 
integer (See Appendix C.) 
integer (See Appendix C.) 
functions (See Appendix C.) 
time 

speed at arbitrary altitude 
speed on entry to the atmosphere 
along range distance 


y 

Y 

Y 0 ( ) 
z 


Kl 

( £ “b) x 

(*■*), 

5 




altitude 

dimensionless altitude, py 

Bessel function of second kind of zero order 
across range distance 
angle of attack 

angle of. attack on entry to the atmosphere 
angle of trim 

density exponential (See eq. (l).) 
half- angle of cone - - . 

deflection target error due to ag 
range target error due to c&g 
deflection target error due to ctp 
range target error due to orp 
altitude variable (See eq. (23b).) 
angle -of- attack function (See eq. (11).) 



MCA TN 1K&8 


29 


angle between flight path and earth’s surface that 
missile has on entering the atmosphere 

air density 

air density at sea level 
radius of gyration 

an arbitrary variable (See Appendix B.) 
oscillation frequency 

Superscripts 

( )* maxim um value of the bracketed parameter which occurs 

in any particular cycle of oscillation 

Subscripts 

( ) maxi mum value of the bracketed parameter 

v 'max * 

Except for the parameters C , f , F , and k: 

( ) x value of the bracketed parameter at altitude for 

maximum deceleration 

( ) value of the bracketed parameter at altitude for 

maximum turbulent heat-transfer rate 

( ) 3 value of the bracketed parameter at altitude for 

maximum laminar heat- transfer rate 

( ) value of the bracketed parameter at altitude for 

maximum normal force 


% 

P 

p o 

a 

f 

Cd 



30 


MCA TN 4048 


APPENDIX B 

DETERMINATION OP ANGULAR VELOCITY AND ANGULAB ACCELERATION 


It is the purpose in this appendix to derive the expressions for the 
angular velocity and acceleration from equation (15b) which gives the 
angular displacement as a function of altitude. Noting that 


dJ 0 (t) 

dy 


- -U(*) f 

a y 


then differentiation of equation (l^b) yields 


d ( a /°®) , „ -py kie 

= -kjJJe 

ay 


-Py 


4 


k 2 + ki e 



P 

2 


hie 


-Pyr 


2 si k 2 +" k x e 


M / f 

J i( 2 4 kg + k x 



(Bl) 


and further noting that 


• q . l>Ji(f)] 


= ^Jo(^) 


dfr 

dy 


then 


a 2 (a/a E ) _ 2 kie' P Y,. 2 . -spy _“Py 


dy e = (ki 2 e~ 2py - k 2 e“ Py j jJ^2 4 k 2 + ki e _ 2 J 


_ JSN 


_ 3Py -py / _ Qf\ 

2p 2 k 1N / k 2 + ki e " 2 "" e^ 1 J x f 2 ^k 2 + e 2 ) 


(B2) 


Using equations (5a), (8a), and (Bl) 



d(a/a E ) 


f. \ 


8v 


dt = PVgBin 9ge v ' ' kie' Py Jo(a^l£e+ ki e“ ^ - 4 k 2 + k x 


and using equations (5b), (8a), (Bl), and (B 2 ) 


By 

2 J 


^2 ^k 2 + kj_ 


py\n 

2 


(B 3 ) 


d^Ctt/qg) 

dt 2 


= P 2 V E 2 ain 2 0 E e (kl ~ ^ - ^)e' 2Py - k 2 e _Py Jjd(Vk 2 + k x e" 


(ski - ~)jkTT 

\ - a / ■ - 


3 By 


k n e 2 J -il 2 J ks> + k 
A ' “ 


^Jkp. + ki e 
A " /J 


(Bin 

\ / 


Since*, for the cases of usual interest, the coefficient k a is so very large compared to k^ 
then one may use the approximate expressions 


_ - fiy (W - / pa 

— - y t ■' - -*Ae BVgSin 0ge 2 e^ / e" 


(B 5 a) 


which, for large values of the argument "becomes 


b ) _ j v/ka 

at - -j n 


BV^sin 0T7.e 


M (kx - ^V Py / _ By\ 

^ e ^ ' cosf4 — 2JkI e ) 




so that idle maximum angular velocity during a cycle is 


u> 

H 
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pw ■ y? ? e( ki ■ ^ py 


Similarly, since k 2 is very large compared to J or kx then approximately 


d 2 fn./n™^ 


\ -By / 


py\ 




which for large values of the argument becomes 


d 2 (a/a E ) ke 4 


p 2 V E 2 sin 2 e E e 


at 2 

% 

and the maximum angular acceleration during a cycle is 


spy . x -By 

“4" ^o)® 


c°b(| - 2s/ki e 




L ct- J ^ JC 


(B5c) 


U) 

ro 


(B6a) 


(Bfib) 


(B6c) 


. 1 , ■ I i ■ I. 
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APPENDIX C 


NUMERICAL INTEGRATION METHOD FOR THE INTEGRAL Ftk^k^ka) 


The following method for the numerical solution of the function 
F(ko,kx,k 2 ) was devised by Dr. William Mersman of Ames Aeronautical 
Laboratory. The integral to be evaluated is , from equation (24a) , 


kaF(ko,k-L,k 2 ) 


-r 




Io(^c)dC 


(Cl) 


where 


and 


a — 2 \j k 2 
b = ko/8k 2 



k^)/4k2 


Io(e,c) 




1 Si J 0 (5 x )46x 


(C2) 


(C3) 


By reference J, page 45, 


5 


r+i 


J r(^) 



r = 0,1,2,. 


(04) 


Introducing 




r+i 




r = 0,1,2,. . . 


(05) 


an integration by parts gives the recursion formula 
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I r (^c) = e _c ^^ r+1 J r+1 (0 + 2c e‘ c ^i 2 ^ ;L r+2 J r+1 (t ;L ) 


That is, 


I r (£,e) = e _c ^ 2 ^ r+1 J r+1 (0 + 2cl r+1 (£,c) , r = 0,1,2 (C6) 


If this equation is multiplied by (2e) r s and then summed on r 
from s to infinity, the following series representation is obtained 


00 

Is(5,e) = e“ C ^^ s+1 £ (2c£) h J 8+h+1 (£) 


(C7) 


h=o 


for any s = 0,1,2,. . ., the series being convergent provided that 

|2c£ | < |. 

In particular, then, setting s = 0 gives 


00 

Io(^c) = e“ c ^ ^ (2cO h J h+1 (?) 


h=o 


and substitution in equation (Cl) gives 


00 

k 2 F(ko,ki,k 2 ) = S s (a,b,c) 


s=o 


(08) 


where 


S s (a,b,c) = (2c ) G f e“^ c_b ^ 2 ^ s J s+1 (C)d^ , ' s =0,1,2,. . . (C9) 


the series is convergent • if 2ca <. 
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Consider, first, tlie term s = 0 


S Q (a,b,c) = f e Ji(?) d 5 

o 

By reference 7 , page 18, Ji(0 = -dJ Q (^)/d^. Hence integration by parts 
gives 

S Q (a,b,c) = J o (0) - e' (c_b)a 2 J 0 (a) - 2(c-b) f e‘ (c “ b) ^J o (0 d £ 


Referring to equation (C3)., this is, since J o (0) = 1, 

S 0 (a,b,c) = 1 - e b ^ a J 0 (a) - 2 (c-b)l 0 (a,c-b) 


(CIO) 


To obtain a recursion formula for the general S B , substitute in 
equation (C9) the equation 


WO -j- J a (0 - 


from reference 7, page 45. This gives Immediately 

S s (a,b ,c) = itcsSg.^a,!^) - (2c) S I s _ 1 (a,c-b) , s = l,2,3,. . . (CU) 

Thus, equations (c8), (CIO) , and (Cll) reduce the problem to one of com- 
puting the sequence I s (a,c-b), s = 0 , 1 , 2 ,. . . . 

For computing purposes , it is desirable to introduce slowly varying 
quantities. The following substitution turns out to be convenient 

P r (a,c-b) = e + ^ C b ^ a I r (a,c-b)/a r+1 r = 0,1,2,. . . (C12) 


The computing problem is then summarized by the following formulae: 
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P s (a,c-b) = Y, [2a(c-b)] h J s+h+1 (a) (C13) 

h=o 

for any s = 0,1,2,. . . . 

P r (a,e-b) = J r +i( a ) + 2a(c-b)P r+1 (a,c-b) , r = 0,1,2,. . . (Cl4) 

S 0 (a,b,c) = 1 - e [J 0 (a) + 2a(c-b)P 0 (a,c-b) ] (C15) 

S s (a,b,c) = 4csS s _ 1 (a,b,c) - (2ca) s e _ ^ c_b ^ a2 P s _ 1 (a,c-b) , s =1,2,3,. . 

(Cl 6) 

00 

k 2 F(k 0 ,k:L,k 2 ) = ^ S s (a,b,c) (C8) 

s=o 


The order In which the computations are to be performed is dictated 
by the following inequalities, each of which is obvious from the corre- 
sponding integral definition, ' under the following general assumptions 


b > 0 , c > 0 , c-b^O 
2ca < 1 , 2a(c-b) < 1 


The inequalities are 


P s (a,c-b)( 


< 


1 

1 - 2a(c-b) 


s = 0,1,2,. 


(017) 


|S B (a,b,c)| (2ac) 8 

1 s + 1 


s =0,1,2 


(C18) 


From the latter, it can be determined how many terms of the series (C8) 
are needed for any desired accuracy, so 'that (C8) is replaced in 
practice by 
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fcfiF(ko,ki,Se) 


K 

^ S s (a,l3,e) 
8=0 


(C19) 


In the present work K = 20, giving a truncation error of less than 10 -10 . 
Once K is chosen, P K is obtained from equation (C13) with s = K where 
again enough terms are taken to insure the desired accuracy. In the pres- 
ent work the computer automatically continued the series (C 13 ) until the 

summand [2a(c-b) ]^Jg- + h +1 (a) became less than 10“ 9 . Once Pg has been 

obtained from the series (C13), the recurrence relations (Cl4) are used 
to compute Pr-i, ?k- 2 > * * • Pi* Po i n that order. Then S D is computed 
from equation (C15), and Si, S 2 ,. . . , in that order from 
equation (Cl6). 

In the present paper the significant range of parameters is 


0 ^ ko ^ 10 

-80 < k x < 0 
5X10 4 < k 2 < 8xl0 5 


For this range it can be shown that the series (C19) can be truncated at 
K = 1 with an error less than lxlO -6 . Furthermore, in the expression for 
S Q and Si, equations (C15) and (Cl6), the terms involving the exponential 
function are also less than lxlO -8 . This gives the simple approximate 
formulae . . . . 

Sq = 1 , Si = 4c 


and, hence 


F(ko,ki,k 2 ) = -i- 
k 2 


1 + 


2 “ 


(C20) 


with an error less than 2xl0 -6 /k 2 . 
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APPENDIX D 

ORDER OF MAGNITUDE OF FACTORS AFFECTING STABILITY 


•In the analysis of this report, a number of simplifying assumptions 
were made regarding the relative importance of the severaljfactors which 
influence the stability. It is the purpose. herein to demonstrate that_ 
three of the assumptions which are of particular importance are, in fact, 
justified.. These. -assumptions are: - 

(1) In the evaluation of values of k 2 (eq. (10c)) 
it is permissible to ignore the Cj_^ contribution 

in comparison with the C^ contribution. (This 
assumption is desired but - not required.) 

(2) In terms involving k 2 + k x , that kj. is 
unimportant. (This assumption is desired but not 
required . ) 

(3) That k 3 - kj 2 is trivial in comparison with 
k 2 + k! (or, from assumption (2) ,^JLn comparison with 
k 2 ) . (This assumption is required to obtain the 
solution for the fundamental differential equation of 
motion (eq. (12).) 

It is to be noted at the outset that k x , k 2 , and k 3 can be written 
in terms of the drag parameter, kQ, in the forms 



The demonstration of the validity of the assumptions (l), (2), and (3) 
will be considered in the sections I , II , and III as follow for conical 
shapes . 


I. For this demonstration it is necessary to show that 
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(3Z, sin 


®E 


is trivial in comparison with 


fWif 
C D W 


Since the Cp^ term is larger the larger the value of l and sin 0 E , 
it will be assumed that the 0 -g is about (slightly greater than) 4 - 5 ° 
and the length is that for a 6 -foot diameter base (i.e., I = ( 3 /tan 5 ) ft). 

Then (since p -1 = 22,000 ft) the comparison of the components is between, 
appr oximat ely 


and 


Cptan 8 



Assuming Z C gA i B 0.30, 0.50, and 0.70 then the ratio of the exact value 
of k 2 to the approximate value of k 2 obtained by ignoring the Cp^ 

term is that shown in figure 16. It is seen that the approximation is 
excellent except when goes to zero (shown for the l C g = O.'JO case). 

Of course, in no practical case would be allowed to approach zero. 


II. To show that k x 
necessary to show that 


Is trivial in comparison with k 2 it is 


ki _ 1 I", Cl a , + 2^ /Vf' 

ko 4 L CD V C D Kv _ 


kg 


1 

2 pl sin 0 g 





is small compared to 
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Again the quantity k 2 will be least compared to k x when l sin 0 E 
is largest; hence, the length and angle assumption of section I is used 
so that the comparison will he between 


ki 

k<3 


and 


1 

k 


kg = 

k 0 



Ass umi ng I cg /z = 0.30, 0.50, and 0.70 then the ratio of the value 
k x to the value k 2 is that shown in figure 17. It is seen again 
Cr 


k £ 

that, except when goes to zero, the approximation is excellent. 

III. In this demonstration it is necessary to show that 


(k 3 


kl 2 )e- 2 Py 



is trivial compared to 


(k 2 + k x )e-£y = k e e-py * kd^e'py 

The test of this assumption is more severe when ko is largest and y 
is minimum, and again, when l sin 0 E is largest. Using the largest 
k Q to be expected, say 20, and the length and angle from section I, and 
for y = 0, the severe comparison Is then between 

*3 ■ * to0 [|% - (|) . 


k 2 


tan 5 /QmqX / 1\ 

10" 5 Ww 


and 
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Assuming l C g/z = 0.30, 0.50, and 0.70, then the ratio of (ks-ki 2 )/^ 
is that shown in figure 18. It is seen again that, except when Cn^ 
goes to zero, the approximation is excellent. Moreover since these coef- 
ficients enter the differential equation (12) multiplied by e" 2 3y and 
e”Py, respectively, and these exponential values have, at altitude, the 
values 


Y, 

ft 

e- 2 Py 


0 

50,000 

100,000 

150.000 

200.000 
250,000 

1.00 

i.o6xio“ 2 
1 . 12 X 10 “' 4 
1.19XIO" 6 
1.25x10“ 8 
1.31x10“ 10 

1.00 

1.03x10“ 1 

1 . 06 X 10 -2 

i.09xio“ 3 

1 . 12 X 10 -4 

1.15x10“ 5 


then the integrated influence of the term involving k 3 - k x 2 
trivial compared to the k 2 term. 


muBt be 
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TABLE I.- ERROR FUNCTION, F x (ko) , EQUATION (26d) 


^o 

Fi(ko) 

0 

0.500 

1 

*571 

2 . 

.659 

3 

• 773 

4 

.921 

5 

1.116 1 ' 

6 

1.376 

7 

1.728 

8 

2.208 

9 

2.873 

10 

3.800 


Note: The values should not "be used when k Q is such that 




TAHLE H.- DYNAMIC STABILITY PARAMETER, t^ko, FOR CONES 




-59-74 

-l4.6l 

-6.256 

-3.339 

-1-997 

-1.276 

-.8510 

-.5856 

-.4168 

-.3102 

-.2476 

-.2183 

-.2149 

-.2317 

-.2617 


73. 00 

iA oA 

8.145 

4.596 

o o=n 

2.051 

1.504 

1.143 

.8890 

.6985 

.5* 18 5 

.4233 

•3147 

m 7rt 
•^17 



0.2 

0.3 

-58.87 

-58.52 

-14.46 

-l4.4l 

-6.243 

-6.254 

-3.376 

~3.4o8 

-2.058 

- 2.102 

-1.353 

-1.402 

-.9380 

-.9930 

-.6818 

-.7386 

-.5188 

-.5770 

-.4157 

-.4730 

-.3530 

-.4070 

-.3188 

-.3672 

-.3050 

-.3450 

-.3046 

-.3347 

-.3122 

-.3313 



-58.57 

-14.46 

- 6.300 

-3.457 

-2.153 

-1.456 

-1.048 

-.7950 

-.6318 

-.5245 

-.4534 

-.4067 

-.3761 

-.3567 

-.3449 



-3.532 

-2.206 

- 1.500 

-1.088 

-.8320 

-.6668 

-.5563 

-.48io 

-.4293 

-.3934 

-.3688 

-.3520 


TABLE HI.- STATIC STABILITY PARAMETER, 


81.00 
on ok 

*-\/ • J • 

9.090 I 10.21 

505© I 5.820 

300 I -3 "7-711 
" J • -J“- ■ ..... 

2.320 
1.708 
1-301 
1.011 


.1328 .1388 



100.7 

110.1 

110.9 

or )iA 

oA i nA 

oA TO 

11.50 

12.82 

13-48 

6.605 

7.445 

7-530 

k 3iR 

1; _ ml; 

R knR 

3-046 

3.492 

3.889 

2.255 

2.591 

2.899 

1.717 

1.967 

2.197 

1.324 

1-505 

1.666 

1.021 

1.145 

I.252 

77W 

Ar£a 

noon 

.5715 

.6195 

.6565 

.4014 

.4264 

.4449 

_p faft 

ryr\ r 

0*7 Ok 

.1487 

.1524 

.1550 


- 1.493 

-1.072 

-.8170 

-.6548 

-.5482 

-.4760 

-.4265 

-.3919 

-.3682 

-.3518 


76.65 

Ol 70 



0-7 

0.8 

0.9 

1.0 

-77.59 

-90.02 

-90.13 

-85.16 

-18.19 

-20.46 

-20.51 

-19.59 

-7.373 

-7.859 

-7.813 

-7.547 

-3.728 

-3.696 

-3.567 

-3.454 

-2.151 

-1.956 

-1.769 

-1,67a 

-1.370 

-1.139 

-•9240 

-.8060 

-.9510 

-.7321 

-.5144 

-.3723 

— • 7120 

-.5238 

-.3202 

-.1654 

-.5702 

- .4167 

-.2380 

-.0833 

-.4830 

-.3631 

-.2142 

-.0708 

-.4277 

-3323 

-.2194 

-.0956 

-.3923 

-.3272 

-.2385 

-.1388 

-.3691 

-.3246 

-.2622 

-.1811 

-.3539 

-.3259 

-.2857 

-.2357 

-.3439 

-.3284 

-.3057 

-.2765 

s( W> 

FOR CONES 


-44.85 

-203.6 

-243.8 

-211.9 

_k RkA 

_kn rA 

■ * w« /v 

= ro on 

_kA no 

1.992 

-H.72 

-18.24 

-17.94 

3.466 

-2.837 

-6.850 

-7.715 

■3 R33 

o£An 

« Jwwu 

1 OR 

_Q Q/V 7 
~*J J 


2.298 

1.746 

1.306 

OROR 

.6755 

.4541 

oA *31 
**- v .JM* 

.1562 


3.148 

2.636 

2.126 

1.666 

I.271 



1.481 

1.8o4 

1.639 

1.403 

1.131 

AAon 

* w ( yu 

.6350 
.4367 
.2 7 68 
.1543 


























Altitude , thousands of feet 
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Cone half-angle , 8 , deg 
(c) Body volume. 

Figure 4.- Continued. 









Dynamic stability parameter 
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Figure 5. -Dynamic stability parameter for example conical 

missiles. 







Altitude, thousands of feet 
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Figure 7. -Altitude variation of oscillations. 




Attitude, thousands of feet 




Oscillation amplitude ratio at altitude for maximum 
laminar heating rate , ( a /a £ )* 
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Cone half- angle, 8, deg 


Figure 8. - Amplitude ratio at altitude for maximum laminar 

heating rate. 


Altitude, thousands of feet 
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Ratio of heat rate at y to heat rate at y 

O 


0 .05 .10 .15 .20 .25 

Oscillation amplitude ratio, ( a /a E )* 

Figure 9.- Laminar heating rate and amplitude ratio as 
a function of altitude for k 0 = 5 , S = 25°. 



Normal force per degree initial angular misalignment 
in terms of missile weight 
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Figure 10.- Maximum ' normal ~ force experienced by example 

missiles. 




Altitude, thousands of feet 
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Figure II.- Altitudes for maximum normal force and drag 

for example missiles. 



Range target error per degree initial angular 
misalignment, feet per deg 
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Cone half -angle, 8, deg 


Figure 12. - Range target error per degree initial angular 
misalignment for example missiles. 



Range target error per degree trim angle, 
millions of feet per deg 
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Figure 13.- Range target error per degree trim angle for 

example missiles. 





Maximum frequency of oscillation, cycles per sec 
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Rim speed for spin rate equal to oscillation 
frequency , feet per sec 
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x 



Cone half- angle, 8 , deg 


Figure 15.- Rim speed for example missiles if spin rate 
equals maximum oscillation frequency. 









Ratio of ( k 2 + k, ) to 


P 
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Figure (7. -Ratio of ( k 2 + k ( ) to k 2 for example of 
largest probable error as a function of missile -cone 
half - angle. 




